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The dirac quasinormal modes of the Reissner-Nordstrom black hole surrounded by quintessence are investi- 
gated using the third WKB approximation. We find that the magnitude of the imaginary part of the quasinormal 
mode frequencies increases firstly and then decreases as the charge Q increases, but it decreases as the absolute 
value of w q increases. The magnitude of the imaginary part of quasinormal complex frequencies is smaller than 
those with no quintessence. That is to say, the dirac field damps more slowly due to the presence of quintessence. 
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' I. INTRODUCTION 

+- > 

Q I The concept of QNMs (quasinormal modes)!! B II & B [1 0, 0, & El El El El El El El El was firstly pointed out by 
Vishveshwara[ 1 8] in calculations of the scattering of gravitational waves by a black hole. Due to the waves emitted by perturbed 
black holes there are no normal mode oscillations but instead of quasi-normal mode frequencies, with the real part representing 
the actual frequency of the oscillation and the imaginary part representing the damping. Quasinormal mode frequencies are 
' the solutions of the perturbation equations, which satisfy the boundary conditions appropriate for purely ingoing waves at the 
horizon and purely outgoing waves at infinity. Quasinormal modes are regarded as the "characteristic sound" of black holes. 
Indeed, quasinormal ringing could provide the direct evidence of the existence of black holes if observed by LIGO (the Laser 
?_i ' Interferometric Gravitational Wave Observatory), GEO600, VIRGO, TAMA300, etc. 

Recently, Kiselev lll9ll considered Einstein's field equations for a black hole surrounded by the quintessential matter and obtained 
a new solution, which depend on the state parameter vv^ of the quintessence. From then on, several people investigate the property 

£SJ , of the black hole suit ounded by quitessenceH Hfl, EE SI H SI M, Ezl 3 , and get some valuable conclusion. In thi s paper , 
we use the third-order WKB approximation to explore the quasinormal modes of massless Dirac perturb ation |29[ [3QL l3~ll l32l [33ll 

OO in RN(Reissner-Nordstrom) black hole surrounded by quintessence. 
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II. DIRAC EQUATION IN THE SPACETIME OF THE REISSNER-NORDSTROM BLACK HOLE SURROUNDED BY 



QUINTESSENCE 

The metric for the RN black hole surrounded by the static spherically-symmetric quintessence can be expressed asJllEl 



ds 2 = -fdt 2 + f- l dr 2 + r 2 (d6 2 + sin 2 6d(f> 2 ), (1) 



with 



2M Q 2 c 

where M represents the black hole mass, w q is the quintessential state parameter, c is the normalization factor related to p q = 
~f zmkp > an d Pq is th e density of quintessence. 

The massless Dirac equation in a general background spacetime can be written as lK34ll 

[Z^ + r^^o, 0) 

where y are the Dirac matrices, 
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e a is the inverse of the tetrad e°, which defined by the metric g^, 



with r\ a b = diag(-l, 1,1,1) being the Minkowski metric, 
is the spin connection which is defined as 



(5) 



where e£ v;jU = d^e^ - T"eb a is the covariant derivative of e\, v , and F^ v is the Christoffel symbol, 
we take the tetrad as 



= diag( -\ff, —— , r, r sin 6) . 

yf 



Then, the Dirac equation ® becomes 



V7 



dt 



r- , / d 1 1 df\ y 2 d 1 



y 3 d*¥ 
r sin 8 dip 



Defining 

Eq. ([8]) becomes 



y = /-?«x>, 



y° <9<D r- , / 3 1\ ^ y 2 3 1 



r sin 3y 



Introducing a tortoise coordinate change 



and the ansatz 



with spinor angular harmonics 



b+ jm 



(D = 



'-It 



(4 



l+^+m ym-1/2 
2/+1 Z l 



\ l+\-m vm +ii2 

y 2/+i I i ' 



(for j = l+ -), 



(6) 



(7) 



(8) 



(9) 



(10) 



(11) 



(12) 



Eq.dTOb can be written in the form 



21-1 ± l 



-1 ± l 



(for j = l- -), 



(13) 



The cases for (+) and (-) in the functions F ± and G ± can be put together, and then the two decoupled equations can be expressed 
in the form 



d 2 F 
~drj 
d 2 G 
~drj 



+ (a/-V l )F = 0, 
+ (w 2 -V 2 )G = 0, 



(14) 
(15) 
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with 



v < = ^^^*y£- { )- [ for k "j*'r imdj = l+\\, (16) 

The potentials V\ and V2 are supersymmetric partners derived from the same superpotential. It has been demonstrated that 
potentials related in this way have the same quasinormal mode complex frequencies ll35ll . We shall concentrate just on Eq. [14] 
with potential V\ in evaluating the quasinormal mode frequencies for the massless Dirac field by the third WKB approximation. 



III. WKB APPROXIMATION AND MASSLESS DIRAC QUASINORMAL MODE FREQUENCIES IN THE 
REISSNER-NORDSTROM BLACK HOLE SURROUNDED BY QUINTESSENCE 



WKB method, a semianalytic technique for determining the quasinormal mode complex frequencies of black holes was 
presented by Schutz and Will[36] at the second order. Later, S. Iyer and C. M. Will[37] developed the method to the third 
order and R. A. Konoplya l38ll extended it to the sixth order. Generally speaking, for the low-lying QNMs, the approach is 
quite accurately. The method has been used to calculate the quasinormal frequencies for the Schwarzschild[39], Reissner- 
Norstrom fl40ll . Kerr 114 ill . Kerr-Newman[42] black holes, etc. The formula for the complex quasinormal frequencies <y(for the 
third order WKB method) is given by 



0? = [Vo + (-2V^') 1/2 A] - i(n + I)(-2^') 1/2 (l + O) 



(18) 



where 



A = 



y 1 



{-iv^yi 2 1 8 







V" 



Q 



-J—{J-(Yl) (77 + 188a 2 ) 
-2V" l 6912 \V" ) 



1 

384 
1 

288 



/' 



2t/(4)\ 



(V"' 2 V { 
V" 3 



/ y (4)^ 



(51 + 100a 2 ) + 



2304 







V" 



(67 + 68a 2 ) 



(V"'V 
V" 2 



(5)\ 



(19 + 28a 2 ) - — 
v ' 288 



t t/(6) \ 



V 







(5 + 4a 2 )} 



(19) 



and 



a = n + — , V, 



(n) 



d"V \ 



(20) 



We take M = 1, c = 0.001 and M — 1, c — for our calculation. And c — means there is no quintessence. 
FigQ] shows the variation of the effective potential for the massless Dirac field with \k\ in RN black hole surrounded by 
quintessence for fixed w q = -0.6, Q - 0.5 and c - 0.001. From FigQ] we can find that the peak value of potential barrier 
gets higher and the location of peak moves along the right as \k\ increases. Figf2] shows the variation of the effective potential 
for the massless Dirac field with Q in RN black hole surrounded by quintessence for |*j = 5, w q = -0.6 and c = 0.001. From 
Fig|2]we can see that as Q increases the peak value of potential barrier gets higher and the location of peak moves along the 
left. Fig|3]shows the variation of the effective potential for the massless Dirac field with w q in the RN black hole surrounded by 
quintessence for \k\ — 5, Q — 0.5 and c = 0.001. From Figf3]we can find that as the absolute value w q increases the peak value 
of potential barrier gets lower. 

The quasinormal frequencies of the massless Dirac particles in the Reissner-Nordstrom black hole surrounded by 
quintessence(for fixed c=0.001) or no quintessence for different Q, w q , \k\ and n (for < n < /c|)are presented in tables ITHVTl and 
figures H]|7]by using the third WKB approximation. 
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CONCLUSION AND DISCUSSION 

We have calculated the low-lying QNM frequencies for massless Dirac field in the Reissner-Nordstrom black hole surrounded 
by quintessence using the third WKB approximation. The results indicate that the quasinormal frequencies not only depend on 
the mass M, the charge Q, but also the quintessential state parameter w„. 
From the tables HITVll and figures|4]|7]we find the main results as follows: 

a)The magnitude of the real part of the quasinormal mode frequencies increases, while the magnitude of the imaginary one 
increases firstly to a maximum at the vicinity of Q = 0.7 and then decreases as the charge Q increases for fixed c - 0.001, \k\, n 
and w q . b)The magnitude of the real and the imaginary part of the quasinormal mode frequencies decreases as the absolute value 
of w q increases for fixed c = 0.001, \k\, n and the charge Q. It means that when the absolute value of w q is larger, the oscillations 
damp more slowly. c)The magnitude of the real and the imaginary part of quasinormal mode frequencies with quintessence are 
smaller than those with no quintessence for fixed |a:| and n. Therefore, we conclude that due to the presence of quintessence, the 
oscillations of the massless Dirac field damps more slowly. 
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FIG. 1: Variation of the effective potential for the massless Dirac field with \k\ in RN black hole surrounded by quintessence for w q = -0.6, 
Q = 0.5 and c = 0.001. 




FIG. 2: Variation of the effective potential for the massless Dirac field with Q in RN black hole surrounded by quintessence for |k| = 5, 



w q = -0.6 and c = 0.001. 
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FIG. 3: Variation of the effective potential for the massless Dirac field with w q in the RN black hole surrounded by quintessence for \k\ = 5, 
Q = 0.5, c = 0.001. 



TABLE I: The quasinormal frequencies of massless Dirac perturbations in the RN black hole without quintessence(c=0) for \k\ = 1, \k\ = 2, 
1*1 = 5 



M Q 



oj(n = 0) 



M Q 



a>(n = 0) 



a>(n = 1) 



0.1 
0.2 
0.3 
0.4 
0.5 
0.6 
0.7 
0.8 
0.9 



0.17682 - 
0.17795 - 
0.17989 - 
0.18274 - 
0.18668 - 
0.19195 - 
0.19895 - 
0.20829 - 
0.22069 - 



0.99 0.23364 



0.10015/ 
0.10028/ 
0. 10048/ 
0.10074/ 
0.10102/ 
0.10125/ 
0.10125/ 
0.10059/ 
0.09790/ 
0.09022/ 



0.1 
0.2 
0.3 
0.4 
0.5 
0.6 
0.7 
0.8 
0.9 



0.37928 - 
0.38128 - 
0.38472 - 
0.38976 - 
0.39668 - 
0.40593 - 
0.41821 - 
0.43479 - 
0.45808 - 



0.99 0.48907 



0.09660/ 
0.09675/ 
0.09701/ 
0.09735/ 
0.09774/ 
0.09812/ 
0.09836/ 
0.09811/ 
0.09632/ 
0.08978/ 



0.35432 - 
0.35649 - 
0.36023 - 
0.36573 - 
0.37330 - 
0.38344 - 
0.39693 - 
0.41503 - 
0.43964 - 
0.46679 - 



0.29888/ 
0.29927/ 
0.29989/ 
0.30070/ 
0.30158/ 
0.30232/ 
0.30247/ 
0.30092/ 
0.29435/ 
0.27401/ 



w(n = 0) 



coin = 1) 



u)(n = 2) 



a>(n = 3) 



ca(n = 4) 



0.1 
0.2 
0.3 
0.4 
0.5 
0.6 
0.7 
0.8 
0.9 
0.99 



0.96183 - 
0.96674 - 
0.97516- 
0.98753 - 
1.00450 - 
1.02716- 
1.05731 - 
1.09783 - 
1.15585 - 
1.23577 - 



0.09631/ 
0.09647/ 
0.09672/ 
0.09707/ 
0.09747/ 
0.09786/ 
0.09813/ 
0.09857/ 
0.09630/ 
0.08994/ 



0.95123 ■ 
0.95620 - 
0.96474 ■ 
0.97727 ■ 
0.99449 ■ 
1.01751 ■ 
1.04814- 
1.09018 ■ 
1.14814 ■ 
1.22708 ■ 



0.29034/ 
0.29080/ 
0.29155/ 
0.29256/ 
0.29372/ 
0.29486/ 
0.29558/ 
0.29649/ 
0.28974/ 
0.27045/ 



0.93165 - 
0.93674 ■ 
0.94548 ■ 
0.95832 ■ 
0.97599 
0.99963 
1.03113 ■ 
1.07571 ■ 
1.13360- 
1.21021 ■ 



0.48789/ 
0.48864/ 
0.48986/ 
0.49146/ 
0.49330/ 
0.49507/ 
0.49608/ 
0.49643/ 
0.48545/ 
0.45273/ 



0.90530 - 
0.91054 ■ 
0.91955 • 
0.93280 ■ 
0.95106 
0.97552 • 
1.00816 
1.05563 - 
1.11354- 
1.18603 - 



0.68959/ 
0.69061/ 
0.69225/ 
0.69441/ 
0.69687/ 
0.69915/ 
0.70027/ 
0.69888/ 
0.68405/ 
0.63751/ 



0.87383 - 
0.87926 - 
0.88860 - 
0.90235 - 
0.92132- 
0.94678 - 
0.98076 - 
1.03093 - 
1.08918 - 
1.15549 - 



0.89485/ 
0.89613/ 
0.89818/ 
0.90087/ 
0.90390/ 
0.90664/ 
0.90777/ 
0.90369/ 
0.88547/ 
0.82508/ 



7 



TABLE II: Quasinormal frequencies of massless Dirac perturbations in the RN black hole surrounded by quintessence for \k\ = 1, \k\ = 2, 



\k\ = 5, and c = 0.001, 



3w q + 1 



-0.01 



W Q 



u>(n = 0) 



W Q 



oj(h = 0) 



a>(n = 1) 



0.1 0.17657 -0.09994/ 0.1 

0.2 0.17770 -0.10007; 0.2 

0.3 0.17963 -0.10027i 0.3 

0.4 0.18248 -0.10053; 0.4 

0.5 0.18640 -0.10081; 0.5 

0.6 0.19165 - 0.10104; 2 0.6 

0.7 0.19864 -0.10104; 0.7 

0.8 0.20795 -0.10039; 0.8 

0.9 0.22030-0.09772; 0.9 



0.37872 - 0.09640; 
0.38071 - 0.09656; 
0.38414 - 0.09681; 
0.38917-0.09715; 
0.39607 - 0.09754; 
0.40529 - 0.09792; 
0.41754 -0.09816; 
0.43406 - 0.09791; 
0.45727 - 0.09614; 



0.99 0.23323 - 0.09010; 0.99 0.48813 - 0.08967; 



0.35381 - 
0.35598 ■ 
0.35971 ■ 
0.36519 ■ 
0.37274 ■ 
0.38285 ■ 
0.39630 - 
0.41434 ■ 
0.43887 ■ 
0.46598 - 



0.29826; 
0.29865; 
0.29928; 
0.30008; 
0.30096; 
0.30170; 
0.30186; 
0.30031; 
0.29380; 
0.27364; 



Q 



oj(n = 0) 



oj(n = 1) 



a>(n = 2) 



a>(n = 3) 



cj(n = 4) 



0.1 
0.2 
0.3 
0.4 
0.5 
0.6 
0.7 
0.8 
0.9 



0.96037 - 
0.96527 - 
0.97367 • 
0.98600 - 
1.00293 - 
1.02553- 
1.05558- 
1.09622- 
1.15379- 



0.99 1.23334 



0.09611; 
0.09627; 
0.09653; 
0.09687; 
0.09727; 
0.09766; 
0.09793; 
0.09774; 
0.09612; 
0.08983; 



0.94980 ■ 
0.95475 • 
0.96327 • 
0.97577 ■ 
0.99295 - 
1.01590 
1.04644 • 
1.08775 • 
1.14609- 
1.22470- 



0.28975; 
0.29021; 
0.29096; 
0.29196; 
0.29312; 
0.29426; 
0.29498; 
0.29430; 
0.28920; 
0.27011; 



0.93027 - 
0.93534 ■ 
0.94406 - 
0.95687 ■ 
0.97449 - 
0.99806 - 
1.02947 ■ 
1.07195 ■ 
1.13158- 
1.20794- 



0.48689; 
0.48764; 
0.48886; 
0.49046; 
0.49230; 
0.49406; 
0.49507; 
0.49361; 
0.48454; 
0.45214; 



0.90398 
0.90920 
0.91819 
0.93141 ■ 
0.94962 - 
0.97401 - 
1.00655 - 
1.05052 
1.11157 
1.18390- 



0.68818; 
0.68920; 
0.69084; 
0.69299; 
0.69544; 
0.69772; 
0.69884; 
0.69633; 
0.68277; 
0.63667; 



0.87259 - 
0.87800 - 
0.88732 - 
0.90103 - 
0.91994 - 
0.94533 - 
0.97921 - 
1.02488 - 
1.08726 - 
1.15354- 



0.89301; 
0.89429; 
0.89634; 
0.89903; 
0.90204; 
0.90478; 
0.90592; 
0.90219; 
0.88382; 
0.82395; 



TABLE III: Quasinormal frequencies of massless Dirac perturbations in the RN black hole surrounded by quintessence for |k| = 1, M 
M = 5, and c = 0.001, 3w„ + 1 = -0.6 



M Q 



a>(n = 0) 



M Q 



ca(n = 0) 



a>(n =1) 



0.1 
0.2 
0.3 
0.4 
0.5 
0.6 
0.7 
0.8 
0.9 



0.17634 - 
0.17747 - 
0.17940 - 
0.18224 - 
0.18616- 
0.19142 - 
0.19840 - 
0.20771 - 
0.22006 - 



0.99 0.23300 



0.09974; 
0.09987; 
0.10007; 
0.10033;' 
0.10062; 
0.10085; 
0.10086; 

0.10022; 

0.09756; 
0.08995; 



0.1 
0.2 
0.3 
0.4 
0.5 
0.6 
0.7 
0.8 
0.9 



0.37820 - 
0.38019 - 
0.38362 - 
0.38865 - 
0.39555 - 
0.40477 - 
0.41702 - 
0.43354 - 
0.45675 - 



0.99 0.48763 



0.09621; 
0.09637; 
0.09662; 
0.09696; 
0.09736; 
0.09774; 
0.09799; 
0.09775; 
0.09598; 
0.08953; 



0.35334 - 
0.35551 - 
0.35924 - 
0.36472 - 
0.37227 - 
0.38237 ■ 
0.39582 - 
0.41385 - 
0.43839 - 
0.46553 - 



0.29767; 
0.29806; 
0.29869; 
0.29950; 
0.30038; 
0.30114; 
0.30131; 
0.29979; 
0.29332; 
0.27320; 



kl Q 



co{n = 0) 



oj(n = 1) 



a>(n = 2) 



u>(n = 3) 



aj(n = 4) 



0.1 
0.2 
0.3 
0.4 
0.5 
0.6 
0.7 
0.8 
0.9 



0.95904 - 
0.96394 - 
0.97234 - 
0.98467 - 
1.00160 - 
1.02420- 
1.05426- 
1.09490- 
1.15248 - 



0.99 1.23207- 



0.09593; 
0.09608; 
0.09634; 
0.09668; 
0.09708; 
0.09749; 
0.09776; 
0.09757; 
0.09596; 
0.08968; 



0.94849 • 
0.95345 ■ 
0.96197 • 
0.97447 ■ 
0.99164 
1.01459- 
1.04513 ■ 
1.08644 ■ 
1.14480- 
1.22344- 



0.28918; 
0.28964; 
0.29040; 
0.29140; 
0.29257; 
0.29372; 
0.29446; 
0.29379; 
0.28873; 
0.26968; 



0.40477 ■ 
0.93408 ■ 
0.94280 - 
0.95560 - 
0.97322 ■ 
0.99679 - 
1.02820- 
1.07068- 
1.13032- 
1.20671 - 



0.09774; 
0.48669; 
0.48791; 
0.48952; 
0.49137; 
0.49315; 
0.49418; 
0.49276; 
0.48374; 
0.45143; 



0.90277 - 
0.90799 - 
0.91698 - 
0.93020 - 
0.94839 - 
0.97278 - 
1.00532- 
1.04929 - 
1.11034- 
1.18273 - 



0.68682; 
0.68784; 
0.68949; 
0.69166; 
0.69412; 
0.69643; 
0.69759; 
0.69513; 
0.68165; 
0.63565; 



0.87144 - 
0.87685 - 
0.88616 - 
0.89987 - 
0.91878 - 
0.94415 - 
0.97802 - 
1.02369 - 
1.08608 - 
1.15243 - 



0.89125; 
0.89254; 
0.89460; 
0.89730; 
0.90034; 
0.90311; 
0.90429; 
0.90063; 
0.88236; 
0.82264; 
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TABLE IV: Quasinormal frequencies of massless Dirac perturbations in the RN black hole surrounded by quintessence for \k\ = 1, \k\ = 2, 



\k\ = 5, and c = 0.001, 



3w q + 1 



-1.2 



W Q 



win = 0) 



M Q 



win = 0) 



win = 1) 



0.1 0.17590 -0.09940/ 0.1 

0.2 0.17702 -0.09953;- 0.2 

0.3 0.17896 -0.09974i 0.3 

0.4 0.18180 -0.10001/ 0.4 

0.5 0.18572 -0.10030/ 0.5 

0.6 0.19098 -0.10055/ 2 0.6 

0.7 0.19797 -0.10057/ 0.7 

0.8 0.20728 -0.09995/ 0.8 

0.9 0.21966-0.09733/ 0.9 



0.37719 -0.09589/ 
0.37919 - 0.09605/ 
0.38262 -0.09631/ 
0.38766 - 0.09666/ 
0.39457 - 0.09706/ 
0.40381 - 0.09746/ 
0.41608 - 0.09772/ 
0.43264 - 0.09750/ 
0.45590 - 0.09576/ 



0.99 0.23265 - 0.08975/ 0.99 0.48685 - 0.08933/ 



0.35245 - 
0.35462 - 
0.35835 ■ 
0.36384 ■ 
0.37139 ■ 
0.38150- 
0.39496 - 
0.41302 ■ 
0.43760 ■ 
0.46483 - 



0.29667/ 
0.29707/ 
0.29771/ 
0.29854/ 
0.29945/ 
0.30025/ 
0.30047/ 
0.29901/ 
0.29262/ 
0.27259/ 



Q 



oj(n = 0) 



oj(n = 1) 



win = 2) 



win = 3) 



win = 4) 



0.1 
0.2 
0.3 
0.4 
0.5 
0.6 
0.7 
0.8 
0.9 



0.95644 - 
0.96135 ■ 
0.96977 - 
0.98213 ■ 
0.99909 - 
1.02174- 
1.05186- 
1.09259- 
1.15030- 



0.99 1.23008 



0.09562/ 
0.09577/ 
0.09603/ 
0.09638/ 
0.09679/ 
0.09720/ 
0.09749/ 
0.09733/ 
0.09573/ 
0.08949/ 



0.94595 ■ 
0.95091 • 
0.95945 ■ 
0.97197 ■ 
0.98918 - 
1.01217 • 
1.04277 ■ 
1.08416 - 
1.14265 - 
1.22149- 



0.28824/ 
0.28870/ 
0.28947/ 
0.29049/ 
0.29168/ 
0.29286/ 
0.29364/ 
0.29304/ 
0.28805/ 
0.26909/ 



0.92657 - 
0.93163 ■ 
0.94037 - 
0.95319- 
0.97084 - 
0.99445 ■ 
1.02590- 
1.06846- 
1.12822- 
1.20481 - 



0.48433/ 
0.48509/ 
0.48633/ 
0.48797/ 
0.48987/ 
0.49171/ 
0.49281/ 
0.49149/ 
0.48261/ 
0.45044/ 



0.90047 - 
0.90567 
0.91467 - 
0.92790 - 
0.94612 - 
0.97054 - 
1.00311 - 
1.04715 - 
1.10831 - 
1.18091 - 



0.68453/ 
0.68557/ 
0.68725/ 
0.68946/ 
0.69199/ 
0.69438/ 
0.69564/ 
0.69332/ 
0.68004/ 
0.63425/ 



0.86928 - 
0.87467 - 
0.88399 - 
0.89771 - 
0.91663 - 
0.94202 - 
0.97592 - 
1.02164 - 
1.08414 - 
1.15071 - 



0.88824/ 
0.88957/ 
0.89167/ 
0.89443/ 
0.89755/ 
0.90043/ 
0.90176/ 
0.89828/ 
0.88027/ 
0.82082/ 



TABLE V: Quasinormal frequencies of massless Dirac perturbations in the RN black hole surrounded by quintessence for \k\ = 1, \k\ 
\k\ = 5, and c = 0.001, 3w„ + 1 = -1.6 



M Q 



win = 0) 



M Q 



win = 0) 



win =1) 



0.1 
0.2 
0.3 
0.4 
0.5 
0.6 
0.7 
0.8 
0.9 



0.17541 - 
0.17653 - 
0.17847 - 
0.18132- 
0.18525 - 
0.19052 - 
0.19752 - 
0.20686 - 
0.21928 - 



0.99 0.23233 



0.09907/ 
0.09921/ 
0.09943/ 
0.09970/ 
0.10001/ 
0.10027/ 
0.10032/ 
0.09973/ 
0.09714/ 
0.08960/ 



0.1 
0.2 
0.3 
0.4 
0.5 
0.6 
0.7 
0.8 
0.9 



0.37604 - 
0.37805 - 
0.38149 - 
0.38655 - 
0.39349 - 
0.40276 - 
0.41508 - 
0.43170 - 
0.45504 - 



0.99 0.48612 



0.09560/ 
0.09577/ 
0.09603/ 
0.09638/ 
0.09680/ 
0.09721/ 
0.09749/ 
0.09729/ 
0.09558/ 
0.08919/ 



0.35149- 
0.35366 - 
0.35740 - 
0.36290 - 
0.37047 - 
0.38060 - 
0.39409 - 
0.41219- 
0.43684 - 
0.46418 - 



0.29569/ 
0.29610/ 
0.29677/ 
0.29762/ 
0.29858/ 
0.29942/ 
0.29971/ 
0.29834/ 
0.29205/ 
0.27212/ 



\k\ Q 



win = 0) 



win = 1) 



win = 2) 



win = 3) 



win = 4) 



0.1 
0.2 
0.3 
0.4 
0.5 
0.6 
0.7 
0.8 
0.9 



0.95347 - 
0.95839 - 
0.96685 - 
0.97926 - 
0.99630 - 
1.01905 - 
1.04929- 
1.09018 - 
1.14811 - 



0.99 1.22819- 



0.09533/ 
0.09549/ 
0.09576/ 
0.09611/ 
0.09653/ 
0.09696/ 
0.09726/ 
0.09712/ 
0.09556/ 
0.08934/ 



0.94306 • 
0.94805 • 
0.95662 ■ 
0.96919 
0.98646 ■ 
1.00955 • 
1.04026 • 
1.08181 
1.14050 
1.21964- 



0.28735/ 
0.28783/ 
0.28862/ 
0.28966/ 
0.29089/ 
0.29212/ 
0.29295/ 
0.29242/ 
0.28752/ 
0.26865/ 



0.92383 ■ 
0.92892 - 
0.93769 ■ 
0.95056 - 
0.96826 
0.99195 ■ 
1.02351 ■ 
1.06620- 
1.12615 ■ 
1.20303 - 



0.48281/ 
0.48360/ 
0.48487/ 
0.48655/ 
0.48851/ 
0.49043/ 
0.49163/ 
0.49042/ 
0.48170/ 
0.44968/ 



0.89791 - 
0.90316 - 
0.91218- 
0.92545 - 
0.94371 - 
0.96820 - 
1.00086 - 
1.04501 - 
1.10635 - 
1.17923 - 



0.68233/ 
0.68339/ 
0.68512/ 
0.68740/ 
0.69001/ 
0.69252/ 
0.69392/ 
0.69178/ 
0.67873/ 
0.63316/ 



0.86694 - 
0.87237 - 
0.88171 - 
0.89545 - 
0.91441 - 
0.93985 - 
0.97382 - 
1.01964 - 
1.08230 - 
1.14914- 



0.88532/ 
0.88667/ 
0.88883/ 
0.89168/ 
0.89491/ 
0.89795/ 
0.89946/ 
0.89623/ 
0.87852/ 
0.81938/ 



TABLE VI: Quasinormal frequencies of massless Dirac perturbations in the RN black hole surrounded by quintessence for \k\ = 1, \k\ = 2, 



\k\ = 5, and c = 0.001, 3w q + 1 = -1.99 



k| Q co(n = 0) \k\ Q oj(h = 0) oj(n = 1) 

0.1 0.17470 -0.09866!' 0.1 0.37434 -0.09527/ 0.35017 -0.29443/ 
0.2 0.17584 -0.09880;- 0.2 0.37637 -0.09543/ 0.35235 -0.29486/ 
0.3 0.17778 -0.09902/ 0.3 0.37984 -0.09571/ 0.35611 -0.29555/ 
0.4 0.18065 -0.09932/ 0.4 0.38494 -0.09608/ 0.36163 -0.29646/ 
0.5 0.18459 -0.09965/ 0.5 0.39193 -0.09651/ 0.36923 -0.29748/ 
1 0.6 0.18988 -0.09993/ 2 0.6 0.40127 -0.09694/ 0.37941 -0.29841/ 
0.7 0.19691 - 0.10002/ 0.7 0.41367 - 0.09724/ 0.39295 - 0.29880/ 
0.8 0.20629 - 0.09947/ 0.8 0.43041 - 0.09708/ 0.41113 - 0.29755/ 
0.9 0.21878 -0.09693/ 0..9 0.45390 - 0.09541/ 0.43589 - 0.29142/ 
0.99 0.23192 - 0.08943/ 0.99 0.48518 - 0.08904/ 0.46340 - 0.27163/ 
Q u)(n = 0) win = 1) a>(n = 2) a>(n = 3) u)(n = 4) 

0.1 0.94899 - 0.09501/ 0.93879 - 0.28635/ 0.91992 - 0.48099/ 0.89442 - 0.67953/ 0.86386 - 0.88142/ 
0.2 0.95396-0.09518/ 0.94382 -0.28684/ 0.92505 -0.48180/ 0.89969 -0.68064/ 0.86930-0.88282/ 
0.3 0.96250 - 0.09545/ 0.95246 - 0.28766/ 0.93387 - 0.48312/ 0.90876 - 0.68244/ 0.87868 - 0.88509/ 
0.4 0.97502-0.09582/ 0.96514 -0.28874/ 0.94683 -0.48488/ 0.92210-0.68483/ 0.89248 -0.88809/ 
0.5 0.99220 - 0.09626/ 0.98254 - 0.29002/ 0.96465 - 0.48693/ 0.94046 - 0.68758/ 0.91151 - 0.89152/ 
0.6 1.01513 - 0.09670/ 1.00580 - 0.29131/ 0.98848 - 0.48896/ 0.96507 - 0.69027/ 0.93705 - 0.89481/ 
0.7 1.04561 -0.09703/ 1.03673 -0.29223/ 1.02023 -0.49031/ 0.99789 -0.69190/ 0.97114-0.89663/ 
0.8 1.08681 -0.09692/ 1.07856 -0.29179/ 1.06318 -0.48928/ 1.04225 -0.69004/ 1.01714-0.89378/ 
0.9 1.14514-0.09539/ 1.13764 -0.28701/ 1.12346-0.48077/ 1.10388 -0.67732/ 1.08005 -0.87655/ 
0.99 1.22576-0.08920/ 1.21727 -0.26822/ 1.20080-0.44894/ 1.17716 -0.63206/ 1.14727 -0.81787/ 
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FIG. 4: The relationship between the real and imaginary parts of quasinormal frequencies for the massless Dirac perturbations in the back- 
ground of the RN black hole surrounded by quintessence for fixed c = 0.001 and no quintessence for c = with the different state parameter 
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FIG. 5: The relationship between the real and imaginary parts of quasinormal frequencies for the massless Dirac perturbations in the back- 
ground of the RN black hole surrounded by quintessence for fixed c = 0.001 and no quintessence for c = with different Q. 
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FIG. 6: The relationship between the real and imaginary parts of quasinormal frequencies for the massless Dirac perturbations in the back- 
ground of the RN black hole surrounded by quintessence for c = 0.001, \k\ = 5, and w q = -0.6 with different n. 
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FIG. 7: The relationship between the real and imaginary parts of quasinormal frequencies for the massless Dirac perturbations in the back- 
ground of the RN black hole surrounded by quintessence (for fixed c = 0.001) and no quintessence (c = 0) for |k| = 5 and Q = 0.5 with 
different n. 



